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We investigate the charge fluctuations of a large quantum dot couoied to a two-dimensional 
electron gas via a quantum point contact following the work of MatveevCtu. We limit our discussion 
to the case where exactly two channels enter the dot and we discuss the role of an anisotropy between 
the transmission coefficients (for these two channels) at the constriction. Experimentally, a channel- 
anisotropy can be introduced applying a relatively weak in-plane magnetic field to the system when 
only one "orbital" channel is open. The magnetic field leads to different transmission amplitudes 
for spin-up and spin-down electrons. In a strong magnetic field the anisotropic two-channel limit 
corresponds to two (spin-polarized) orbital channels entering the dot. The physics of the charge 
fluctuations can be captured using a mapping on the channel-anisotropic two-channel Kondo model. 
For the case of weak reflection at the point contact this has already briefly been stressed by one of us 
in PRB 64, 161302R (2001). This mapping is also appropriate to discuss the conductance behavior 
of a two-contact set-up in strong magnetic field. Here, we elaborate on this approach and also 
discuss an alternative solution using a mapping on a channel-isotropic Kondo model. In addition we 
consider the limit of weak transmission. We show that the Coulomb-staircase behavior of the charge 
in the dot as a function of the gate voltage, is already smeared out by a small channel-anisotropy 
both in the weak- and strong transmission limits. 



I. INTRODUCTION 

In the past few years a great amount of work has 
been deitoted to studying the Kondo effect in mesoscopic 
systemsEl. A motivation for these efforts was the recent 
experimental observation of the^Kondo effect in tunnehng 
through a small quantum dotou. In these experiments 
the effective (or excess) electronic spin of the dot acts as 
a magnetic impurity. 

A different set of problems relating the Kondo effect 
to the physics of quantum dots is encountered when 
studying fluctuations of the charge of a large Coulomb- 
blockaded quantum dot. The setup we have in mind con- 
sists of a large quantum dot coupled to a reservoir via a 
quantum point contact (QPC) and capacitively coupled 
(with a capacitance Cgd) to a back-gate (see Fig. |l|). The 
amount of charge on the dot can be changed through 
the gate voltage Vg- The term "large" implies that the 
spacing A of the energy levels on the dot (almost) van- 
ishes and is much smaller than the dot's charging energy 
Eg = e^/(2Cgd). We assume that the capacitance be- 
tween the two-dimensional electron gas (2DEG) and the 
dot can be neglected. 

There are two limits in which a mapping to a Kondo 
Hamiltonian can be used to calculate the charge on the 
quantum dot. These two limits roughly correspond to 
the cases of very strong and very weak reflection at the 
QPC. Before elaborating on this point we will here give 
the standard multi-channel Kondo Hamiltonian 

Hk = HK^n + J2 [Jo..zS^:MSl (1) 
a 

+ Ja,± {saAO)S'x + Sa,ymSl)] 



in a very general form to be used as a point of refer- 
ence for later discussionsQ. Here S is an impurity spin 
and Sq(0) is the spin of conduction electrons with fla- 
vor a at the place of the impurity. We will consider 
only the two-channel case a = 1,2. In the discussion of 
the Coulomb blockade problem we will encounter Kondo 
Hamiltonians which are both spin- {Ja,z 7^ Ja,±) and 
channel-anisotropic (Ji._l 7^ J2.±)- 




FIG. 1. A large quantum dot is coupled to a 2DEG via a 
quantum point contact. The number of electrons in the dot 
can be controlled through the gate voltage Vg- The auxil- 
iary voltage Va can be used to open or close the point contact 
and thus to adjust the reflection amplitudes for the transport 
channels in the QPC. 

We now want to return to the original problem and 
first want to discuss the limit of weak tunneling between 
the dot and the reservoir (strong reflection at the QPC). 
For low enough temperatures (T <C Eg) the charge on 
the dot is quantized. When the gate voltage Vg is in- 
creased the charge on the dot changes in a step-like man- 
ner. This behavior is referred to as a Coulomb staircase. 
Charge fluctuations are important only for those values of 
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the gate voltage at which two neighboring charge states 
(e.g. with n and n + 1 electrons on the dot) are energy 
degenerate. Matveev demonstrated the equivalence of 
the effective Hamiltonian describing the charge dynam- 
ics close to such a degeneracy point to a spin-anisotropic 
Kondo Hamiltoniand. If the dot and the reservoir are 
connected by m transport channels with the same trans- 
mission probability the problem can be_mapped on a 
channel-isotropic m-channel Kondo modefl with J^^z = 
and Ja,± IQ. The similarity of the charge dynamics on 
the dot to the spin dynamics in a Kondo model was al- 
ready observed by Glazman and Matveev in there discus- 
sion on the conductance of a metallic, grain (a quantum 
dot) tunnel-coupled to two reservoirsQ. 

The limit of a point contact with one (electrons with- 
out spin) or two (electrons with spin) transport channels 
close to perfect transmission has also been treated. An ef- 
fectively one-dimensional model for the dot-QPC system, 
which makes it possible to calculate the total charge o 
the dot using bQSonization,pWas developed by Flensber^ 
and by MatveevQ. MatveevQ also calculated the dot's ca- 
pacitance as a function of the gate voltage and the reflec- 
tion amplitude at the point contact. In his calculation he 
demonstrated that the Hamiltonian of the original prob- 
lem can be mapped onto a Kondo Hamiltonian in the 
generalized Toulouse limit (Ji.z = J2,z = invp)- In this 
limit the two-channel Kondo model is_fixactly solvable 
as was shown by Emery and Kivelsontj. It was one of 
the main results of Ref. ^ that Coulomb blockade os- 
cillations (accompanied by a logarithmic singularity in 
the expressions for the charge or the capacitance close 
to "half-integer" values of the number of electrons in the 
dot) persist also in the limit of weak reflectioii^ as re- 
cently checked experimentally by Herman et aLliil. They 
are completely smeared out only if at least one transport 
channel is at perfect transmission, the charaej.iujbhe dot 
then increases linearly with the gate voltageoElllS. A re- 
view of these results and more generally of interaction 
effects in quantum dots can be found in Ref. |l3| . 

Extending our previous workllil we here want to dis- 
cuss the effect of a channel- anisotropy on the results for 
the capacitance in the two limits of almost perfect trans- 
mission and almost total reflection discussed above. We 
will show that the Coulomb staircase behavior is smeared 
out by a small channel anisotropy both in the weak and 
strong reflection limit. We will concentrate on the case of 
two transport channels. Channel-anisotropy then means 
that the reflection (or transmission) amplitudes for the 
two channels are different. Such a situation can be re- 
alized applying a magnetic field to the sample. A weak 
in-plane magnetic field leads to different reflection am- 
plitudes for spin-up and spin-down electrons. In a strong 
magnetic field electrons are spin polarized and the num- 
ber of open channels in the QPC and their respective re- 
flection amplitudes depend on the opening of the QPC. 
The confinement potential at the QPC and thus its open- 
ing can be tuned changing the auxiliary gate voltage Va 
(see Fig. p. In the limit of small reflection an intu- 



itive physical picture for the transition from the channel 
symmetric caseQ to the case of a weat^asymmetry be- 
tween the two channels can be obtainedEj from the use of 
the channel- anisotropic two-channel Kondo model at the 
Emery-Kivelson linolj. This model has also been investi- 
gated in Ref. |l6|. The two coupling constants of the asym- 
metric Kondo model are directly related to the reflection 
amplitudes of the two transport channels in the QPC. In 
this paper we will elaborate on this mapping and discuss 
its limitations. Although the physical picture associated 
with the channel- anisotropic two-channel Kondo model 
is very appealing, the Coulomb blockade problem in the 
high transparency limit with two transport channels with 
different reflection coefficients can also be solved through 
a mapping on the channel-isotropic two-channel Kondo 
model as we will shoW|_Such a procedure was applied by 
Matveev and FurusakiEJ when they calculated the inelas- 
tic cotunneling through a quantum dot strongly coupled 
to two quantum point contacts with one transport chan- 
nel each (electrons in a strong magnetic field). We will 
see furthermore that in the case of tunnel-coupling (small 
transmission) an asymmetry between transmission am- 
plitudes of different channels will lead to a mapping of 
the Hamiltonian on a channel-anisotropic Kondo model 
with {Ji,z = J2.Z = 0). 

Our paper is structured as follows: In a first part we 
consider the weak coupling limit. We start by develop- 



III we dis- 



ing the model Hamiltonian in Sec. ||. In Sec. 
cuss the mapping on the channel-anisotropic two-channel 
Kondo model and use it to calculate the capacitance of 
the quantum dot in a weak magnetic field. In the fol- 
lowing section (Sec. IV) we show that the same results 
can also be found from a mapping to a channel-isotropic 
Kondo model and in addition look at the case of a strong 
magnetic field applied to the QPC-dot system. In Sec. ^ 
a simple scaling argument is used to rederive the main 
(exact) results of Sees. Ill and IV, In Sec. VI, we will 
show that the channel-anisotropic two-channel Kondo 
model can also be used to discuss the conductance in a 
two-contact set-up with a strong in-plane magnetic field. 
The limit of small transmission through the quantum dot 



will finally be considered in Sec. VH 



II. THE MODEL 



We here consider a large quantum dot weakly coupled 
to a 2DEG via a point contactti (see Fig. |^). The shape of 
the QPC is defined trough metallic gates put on top of a 
2DEG. The lateral confinement potential in the QPC can 
be controlled changing the voltage applied to these gates. 
As a consequence of the lateral confinement the conduc- 
tance of the QPC is quantized and electron motion in the 
vicinity of the QPC is essentially one-dimensional. The 
one-dimensional wave-function ^'JJ for motion along the 
X axis is characterized by the spin index cr =t, i and the 
(orbital) channel quantum number n due to the lateral 
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confinementlla. Note, that from here on we wih use the 
word channel to mean electrons with a certain pair of 
indices cr, n. As a short notation we introduce the chan- 
nel index a = a,n and we denote the wave-function of 
an electron in this channel by '^a- In the further anal- 
ysis we neglect transport channels that are totally re- 
flected at the QPC, since electrons from these channels 
are confined to the reservoir and will not contribute to 
the charging of the dot. We will concentrate on the case 
where only two channels are open. This can be real- 
ized in two different ways. If a weak (or zero) in-plane 
magnetic field is applied (see Sec. Ill), the two channels 
correspond to the two spin polarizations of electrons in 
the lowest (ti = 1) energy eigenstate of the lateral Hamil- 
tonian (^'i = *|,*2 = *|)- The role of the magnetic 
field is to introduce an asymmetry between the reflection 
coefficients for spin-up and spin-down electronscJ. In a 
strong magnetic field the electrons are spin-polarized and 
the two-channel case corresponds to two orbital channels 



To manipulate the Hamiltonian of the one-dimensional 
interacting system we have introduced above, it is con- 
venient to use the bosonization technique. Bosonizing 
in the standard way we express the Fermi field operatocs. 
through the bosonic fields (f>a (x) and 9a (x) and writa23'Eil 



^a,B/L{x) 



1 



(6) 



Here p = I for right- movers (R) and p = —1 for left- 
movers (L). The bosonic fields obey the commutation 
relations 



[(j)a{x),df}iy)] ^ -sgn{x ~ y)Sa,i3, 
[(j)a{x),dyei3(y)] ^iS{x-y)Sa.p, 
[Mx),My)] = [9a{x),9p{y)]=0. 



(7) 



ente ring the constriction (e.g. 
Sec. fv. 



*i = *l*2 - ^t); See 



In the new variables the kinetic energy takes the form 



Since we are interested in energies much smaller than 
the Fermi energy we linearize the spectrum around the 
Fermi points. Furthermore it is convenient to decompose 
the wave-function into a right-going and a left-going con- 
tribution = exTp{ikpx)'^a.R{x)+exTp{—ikpx)'9a.Lix). 
For two open transport channels the Hamiltonian for the 
electrons in the point contact is 



a=l,2 



dx 



(dx4>af + {dxOaf 



(8) 



To bosonizc the Coulomb Hamiltonian we use the rela- 
tion : I'^^'a.L + ^'Lfl^'a.i?, := -a^^a/V^. Thc total 
charge in the dot then is 



HKin — —ivp 



E 



r dx k ^i,)d^^aA-) (2) Q = E I (:■ + 



(9) 



of 



where the sum is over (open) quantum channels. Electron 
interactions in the quantum dot are taken into account 
via the Coulomb Hamiltonian 



He - EciQ - N)^ - EcN^ 



(3) 



The parameter eA^ — VaCgd is proportional to the gate 
voltage and Ec = e^/{2Cgd) is the chargijBg energy. In 
addition we allow for a weak backscatteringE3 in the point 
contact: 



The charge in the dot is now measured in piit 
We neglect both finite size effects in the dotll^E^ (sup- 
posed to be ideal with no dephasing prooeasXand meso- 
scopic corrections to the capacitance Cgd""- Express- 
ing the charge in terms of the bosonic variables we obtain 
Q — (01 (0) + cf)2{^)) /v^- The term at spatial infinity is 
independent of the gate voltage; We choose (j)i{oo) = 
in such a way that the total charge Q on the dot is zero 
when N = 0. Consequently we get 



I <«(0)*„,l(0)-H/j.c. 



Q=l,2 



(4) 



Ho = 



(10) 



The transmission is supposed to be globally adiabaticll3. 
The reflection amplitudes |ri| and \r2\ can be tuned ap- 
plying a voltage to the gates defining the QPC. The main 
goal in the following will be to calculate the shift in the 
ground state energy Se due to the backscattering in the 
point contact. From the correction Se it is possible to 
obtain the average charge in the dot and to dot's capac- 
itance via 



for the Coulomb Hamiltonian. In writing Eq. din) we 
have only considered the Q-dependent part of Eq. (g^. Fi- 
nally we also have to bosonize the backscattering Hamil- 
tonian which leads us to 



Hts - — E 



r„ cos 



(V4^<^„(0) 



(11) 



(0) = CgdVa 



C = 



d{Q) 
dVn 



djSe) 

= Cgd - 



(5) 



dVa^ ' 



Now we introduce the new standard variables (pc s = 
^ {^i{x) ± </.2(x)) and 0,,, = -ij {9i{x) ± 02 (x)) where 

the positive sign belongs to the label c (charge) and the 
negative sign belongs to s (spin) . Note that in the case of 
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a strong magnetic field electrons are spin-polarized. The 
index s then labels a "pseudospin" and not physical spin 
states. The kinetic energy in the new variables has the 
standard form of Eq. (||) . The Coulomb Hamiltonian 



2Ec 



He = ^ (mo) - v^^y 



(12) 



is a function of the charge field only which is the main 
motivation for the introduction of this new set of vari- 
ables. Furthermore we get for the backscattering part 



Hb, = ^(Inl + IrsDcos (^2^0,(0)) cos (V2^0,(O)) (13) 
(In I - \r2\) sin (V2^</.,(0)) sin (^2^0,(0) 



na 

VF , 



We now introduce the charge fluctuation field 0c (0) — 
0c(O) — \PkJ2'N . The total charge in the dot is pinned at 
its classical value Qd = ^f'c(O) = jlN to minimize the 
Coulomb energy. For weak backscattering |ri|,|r2| ^ 1 
and for energies below the charging energy we can aver- 
age over the charge fluctuations 0c (0). When averaging 
over the term cos (A/27r(/)c(0)) — cos ^\/27r0c(O) -t- ttA^ 
in Eq. (RSh we obtain the expression 



-'r(0c(o)'> cos(7r7V) 



Vp 



cos(7rA^). 



(14) 



Here the angular brackets mean averaging with regard 
to the ground-state of the free Hamiltonian Eq. (H). An 
analogous expression can also be found for the term 
sin (\/27r(/)c(0)) in Eq. (p3|). When calculating the cor- 
relator 



:(0) ) = -^^-In 

ZTT Vp 



(15) 



we have taken into account that only fluctuation modes 
with energies larger than the charging energy Ec can 
enter into the dot. Here 7 is defined through ^ = <f 
where C « 0.577 is Euler's constant. After the averaging 
the backscattering part of the Hamiltonian can thus be 
written as 



Eq. (0) will be useful as a starting point for the mapping 
on the isotropic form of the two-channel Kondo model 
(see Sec. pv| ). The charge part of the kinetic energy can 
now be dropped since it is completely decoupled from the 
perturbation due to the backscattering. 



III. CHANNEL-ANISOTROPIC TWO-CHANNEL 
KONDO MODEL 

In this section we will concentrate exclusively on the 
case where the reflection amplitudes |ri| and \t2 \ are very 
close to each other. It is then convenient to introduce the 
parameters = Ij^iI + \t2\ and |i5r| = |r2| — where 
\8r\ ^ \R\ ^ 1. For the sake of clarity we have split up 
this section into three subsections. In the first subsection 
we will discuss the relation between our problem and the 
channel-anisotropic (two-channel) Kondo modelij. In 
the second subsection we will derive an expression for 
the impurity correction to the ground state energy while 
a physical realization of the case \8r\ ^ ^ 1 is dis- 
cussed in the third subsection. There we will also explic- 
itly calculate the charge in the dot and the dot's capaci- 
tance. 



A. Equivalence to the Kondo Model 

We will now see that our theory for the Coulomb 
blockade problem can be mapped on the anisotropic 
two-channel Kondo model at the Emery-Kivelson line 
i.e. t Xi M Y 2TrvF with a = 1,2. The solution of-this 
modelEj'Ea is due to Fabrizio, Gogolin and NozierestH 

Introducing the auxiliary impurity-spin operators Sx 
and Sy we rewrite the backscattering term (Eq. (|l6|)) as 

Hbs = — cos fV2^0,(O)) 5, + ^ sin fV2^0,(O)) Sy, 

(18) 



where the Kondo coupling parameters are defined 
through 



H, 



^/^aEcVp 
Tra 

^/^aEcvp 



(|ri| + Iral) cos {t:N) cos (V2^0,(O)) (16) 
(|ri|- |r2|)sin(7rAr) sin (V2^0,(O)) . 



Jx = 2\Rya'^Ecvp cos(7r7V), 
Jy — 2\5r\ \/ a'yEcVp sin(7riV). 



(19) 



Using simple trigonometric relations the Haniil|tonian 
Eq. ( |l6| ) can be rewritten in an alternative formLj as 



Tra V / 



where r is the complex parameter r = (|ri|e*'^^ + 
|r2|e~*'^^)/2. The first form of the backscattering term 
Eq. (^) will be used in Sec. HI when discussing the two- 
channel anisotropic version of the Kondo model, while 



Starting with a standard Kondo model, we rather get 
Jx oc (Ji,_L + J2.x) and Jy cx {Ji^±_ - J2,_l): where J^-a 
denotes the transverse Kondo coupling of each conduc- 
tion band channel with the magnetic impurity. The to- 
tal Hamiltonian which we will denote H^p, is given by 
H^j^ — HKin{4's, Os) + Hbs{4's)^ where the kinetic term is 
of the form Eq. (||). Both Sx and Sy can be considered as 
good quantum numbers since their commutators with the 
Hamiltonian H^j^ are small close to perfect transmission 
through the quantum point contact: 
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H 



cx 



-i\R\S,, 



(20) 



H^j^,Sy (xi\Sr\S^ 



The impurity spin in H^^^ can oscillate between the 
two values Sx = 1/2 and Sy = —1/2. It is impor- 
tant to note that the backscattering part of the total 
Hamiltonian (Eq. (p^) and the coupling term of the two- 
channel anisotropic Kondo model (Eq. (p^)) are exactly 
equivalent only in the channel symmetric case \Sr\ = 0, 
since Sx and Sy do not commute. However, we will see 
in the following section that the approximation made 
in writing Eq. ( [l^ ) is good as long as \Sr\ ^ \R\ or 
7V-0,±1/2,±1.... 

We now want to find the shift of the energy due to the 
backscattering or, in the language of the Kondo prob- 
lem, the impurity correction to the ground state energy. 
To make further progress it is useful to refermionize the 
Hamiltonian. The basic idea is to introduce a unique op- 
erator tp{x) such that cos(-\/27r</'s(0)) oc 'tp{0) +'^|J^0) and 
sm{^/2n(j)s{0)) oc V'(O) — "0^(0)- Furthermore we use the 
Majorana representation 



(21) 



V2Sx^a = {d + d^)/V2, 
V2Sy ^ -b ^ {d^ - d)/{iV2), 



to express the spin-operators through fermionic opera- 
tors. The d operators {d — a + ib) obey {d^ , d} = 1 and 
{d,ip{x)} = 0. The (unusual) refermionization procedure 
will be extensively discussed in Appendix There we 
also give a precise definition of the new fermionic fields 
^p. After refermionization the kinetic energy takes the 
standard form 



HKin — —ivp 



dx ip\x)dxtp{x), 



while the backscattering Hamiltonian is 
iJx 



i^b. = ^=^(V^(0)+^^(0))& 




(22) 



(23) 



(V;(0)-^t(o))a. 



So we finally arrive at a (standard) solvable resonantJeiiel 
type of model for the two Majorana fermions a and t£j't3. 



B. Impurity corrections and scattering phase shifts 

Our next step is to calculate the impurity corrections 
to the ground state energy from the impurity Green's 
functions Ga = - {T-ra{T)a{{))) and Gb = -(T^6(r)&(0)). 
The Fourier transforms of these Green's functions can be 
found in a convenient way from the equations of motion 
(cf . Appendix P) . We get 



Gk{uj) = 



1 



(24) 



where the label is fc = a, 6. The width of the reso- 
nant a{b) level ra(rf,) is related to the respective Kondo 
coupling constant Jy{Jx) via 



2 



|(5r|^sin^(7rAf), 



(25) 



— 'R\ cos (ttN). 



Airavp TT 

The impurity correction to the ground state energy at 
zero temperature is 



5e 



Ec 



doj 
2^ 



ui (na(w) -I- n}y{uj)) 



(26) 



The occurrence of a high energy cut-off at Ec in Eq. ( [2q ) 
is an intrinsic property of the theory developed so far. 
The low energy cut-off uJmin needs some additional ex- 
planation. We can distinguish two different situations: 
In the neighborhood of integer values of A'' the resonance 
width Ta is negligibly small. Already at temperatures of 
the order of Ft spin fluctuations will get pinned, since 
the coupling constant Jx goes to strong coupling (com- 
pare Sec. 0). 




0.2 0.4 0.6 0.8 



FIG. 2. The resonance energies Va and Vb of Eq. ( |25| ) 
are shown. The parameter |_R| in Tb (full line) is |_R| = 0.4. 
The size of the resonance Va grows with increasing anisotropy 
and thus with increasing \5r\. The lowest (dotted) line in 
the graph corresponds to \5r\ — 0.1 the second lowest (short 
dashes) to \5r\ = 0.2 and the third (long dashes) to \5r\ — 0.3. 
The energies are measured in units of Ec^/ir. 

Very close to half-integer values of N however, we 
have Fb < F^ even though \6r\ < \R\. Thus the 
coupling Jy will go to strong coupling before Jx and 
spin fluctuations will be frozen at Fq. An expression 
for Umin that correctly reproduces these two limits is 
i^min — ?7^ax{Fa, Ffc}. The density of states nk{uj) of 
impurity k in Eq. (|2^) is related to the corresponding 
impurity Green's function through 
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"-fc(^) = -2Ini{Gfc(tj)} = 



2r/cSgn(w) 



(27) 



With this expression for nk{Lij) the integration in Eq. (|2; 
can easily be performed and we get 



(5e 



271 ^ 



max{ra,rb}2 + r;; 



(28) 



For most purposes it is sufficient to. approximate this 
pression by the more simple formEj 



Se 



{Ta + n) In 



Ec 



max{ra,rb} 



(29) 



We have used Tk ^ Eq and have dropped an additive 
constant. A quantity which is interesting is the scattering 
phase shift (extracted from the Friedel sum rule) 



electrons with spin-up and spin-down are different due 
to the Zeeman effect. The channel index a here dis- 
tinguishes between spin-up and spin-down electrons in 
the lowest orbital channel n = 1. The (original) elec- 
tronic wave-functions for electrons in the two channels 
are ^a=i = and *q=2 = *x"^- We first want to 

consider the special case of zero magnetic field which was 
solved by Matveev in Ref. The reflection amplitudes 
for the two channels, corresponding to spin-up and spin- 




R\/2) and thus 
25|) that Jy = 



down electrons are equal (|ri 
|(5r| = 0. It follows from Eqs 
and Ta = 0. Furthermore from Eq. we know that 
Ff, ^\R\'^"fEcCQf'^{TrN)/7T. Calculating the energy with 
Eq. m) we fin. 



cos2(^iV)ln (7/7r|i?|2 cos^ttN)) . (31) 



SiLo) 



1 



y arctan 
2 ^ \ uj 

k—a,b 



(30) 



of the conduction electrons "0 due to the impurity scat- 
tering. In the channel symmetric case we have Fa = 
and 5{u} <^ F^) = 7r/4. The Friedel sum rule can also be 
rewritten Z = 2 {21 + 1) J/tt where Z is the impurity 
charge screened by the electrons. For s-wave scattering 
[l — 0) and Z = 1/2 we recover 5 = 7r/4. The particular 
value of 5 can thus be understood as a consequence of the 
fact that only "half" of the fermion d (the part a) is cou- 
pled to the conduction electrons. In general, if both Fq 
and Ffc are finite we find 5 = 7r/2 in the limit w — > and 
Z = \ since now a and b are screened. Here, (even) for 
a finite channel asymmetry, we still find 5 = 7r/4 at the 
fixed point because close to = 1/2 wc have F;, ^ and 
close to = 0, ±1 we have F^ 0. This gives a physical 
justification on the validity of the mapping in Sec. IV. 
It is important to remember that the "spin-fermions" ij: 
are not the real electrons but are related to the spin de- 
grees of freedom of the original electronic wave- functions. 
However, they play the part of the conduction electrons 
of the real Kondo problem. In the case of spin-less elec- 
trons (electrons in a strong magnetic field) and for a sin- 
gle transmitted channel |r| <C 1 the scattering phase shift 
of the real electrons is related to the average charge on 
the dot via 5 = tt{Q), again as a consequence, of Friedel's 
rule. It was shown by Aleiner and Glazmantj that in the 
one-channel limit this relation can be used to calculate 
the impurity correction (see Eq. (|35|)) in an intuitive way. 
The physics respjibles closely the one of the one-channel 
Kondo problemEII. 



C. Applications 

The case \5r\ ^ \R\ ^ 1 can be realized in our 
setup by applying a weak in-plane magnetic field. In 
a non-zero field the reflection amplitudes |ri| and 17-21 for 



The correction to the capacitance can then easily be cal- 
culated using Eq. (h) and is found to be 



SC = ~2-iEc\RVl3^ cos(27rAf)ln 



|i?|2cos2(7rAr) 



(32) 



Here we have kept only the logarithmically divergent con- 
tribution which will dominate all other terms close to 
A^ = 1/2. The parameter /3 = e/{2Ec) is the ratio of 
the dimensionless parameter A^ and the gate voltage Vg, 
(3 = N/Vg- This seems to be in agreement with the 
recent capacitance experiment of Ref. ^l]. 




FIG. 3. The correction to the capacitance Eq. (^) 
is shown for different values of the anisotropy parameter 
\5r\ — \r2 \ — \ri\. The parameter \R\ = Ij-iI + \r2\ is set to 
|-R| = 0.4. The full line corresponds to the channel isotropic 
case \5r\ — where the capacitance is logarithmically diver- 
gent. In order of decreasing peak height the remaining three 
curves correspond to \Sr\ = 0.1, |5r[ = 0.2 and |5r| = 0.3. 
The units of the capacitance are arbitrary. 
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As we can see from Eq. (jl^) in the channel-symmetric 
case only fermion b is coupled to the field ip. If we al- 
low for a weak magnetic field the reflection coefficients 
for spin-up and spin-down electrons are slightly different 
and |(5r| ^ 0. As a consequence the Kondo coupling Jy 
does not vanish anymore and both Majorana fermions a 
and b are coupled to the bath. The total energy shift 6e 
from Eqs. (|5) and {M) is found to h<M 



5e = +^ cos2(^iV) -I- |5rp sm^inN)) (33) 
X In (7/7rmax{|(5rpsin2(7rA^),|i?pcos2(7rA^)}) . 

Due to the appearance of the max the logarithm does 
not diverge anymore. Still it can be very large due to the 
smallness of |<5r| and the logarithmic term will dominate 
in the capacitance (See Fig. 3): 

dC = -2"fEcl3^ {\Rf - \Sr\^) cos(27riV) (34) 
X In 



max{|5r|2sin^(7riV), |i?|2 cos2(7rA^)} 



We see from Eqs. ( p3| ) and (|3J) that an arbitrary weak 
anisotropy between the two reflection coefficients is suffi- 
cient to cut off the logarithmic divergence. At the degen- 
eracy point N — 1/2, we find SC oc In \Sr\. There is some 
analogyfl with the behavior of the magnetic susceptibility 
of the impurity x = d^{6e)/dh^ away from the Emery- 
Kivelson linal3 and with that of the local magnetic sus- 
ceptibility xi = d{Sz)/dh at the Emery-Kivelson Ipe 
(the magnetic field h would only act on the impurity )E3, 
even though for weak backscattering at the QPC there is 
no real correspondence between the charge in the dot Q 
and Sz- From the two-channel anisotropic Kondo model 
we see that the appearance of the second energy scale 
which leads to the suppression of the divergence follows 
in a natural manner from the coupling of the second Ma- 
jorana fermion to the conduction electrons. 

At this point it is also interesting to compare these 
results to the result obtained in the case of a single 
transmitted channel (reflection amplitude |ri| <C 1 and 
— *■ !)■ SuGh a situation can be realized in a strong 
magnetic fieldllj (spinless electrons) and was also treated 
in Ref. ^ The energy-shift for this case is 



Seoc |ri|£'c cos(27rA^) 



(35) 



and therefore no logarithmic contribution occurs any- 
more in the expressions for the charge and the capaci- 
tance: We only get a periodic oscillation as a function of 
the gate voltage. 

In the next section we want to demonstrate that the 
approximation made in writing Eq. ( p^ is justified. Fur- 
thermore we will address the case of strong asymmetry 
between the conduction channels. 



IV. MAPPING TO A CHANNEL-ISOTROPIC 
KONDO MODEL 

We have seen in the last section that the backscatter- 



ing term Eq. (18) is not exactly equivalent to the orig- 



inal Hamiltonian Eq. (16). We will now justify the ap- 
proximation made and will rederive the results we have 
found from the two-channel anisotropic Kondo model in 
an exact-Jjvay, following in our derivation Furusaki and 
Matveevll3. The Section is structured in a similar way as 
Sec. {fn^. The mathematical mapping will be discussed 
in a first Subsection while a second short Subsection will 
be devoted to the discussion of an application. 



A. Mapping 

As a starting point we use the form Eq. ( p^ of the 
original backscattering Hamiltonian which we rewrite as 



Jo 
2'Ka 



^0^,(0) 



(36) 



Here r is the complex parameter 
^)/2 and Jo 



r = (|ri|e"^ -|- 

|f2|e j/'Z and Jq = A^/ajEcVp- Exactly as in the 
previous section we have introduced an auxiliary im- 
purity spin Sx- The spin Sx is a good quantum num- 
ber since it commutes with the total Hamiltonian (here 
Sx = 1/2). It is important to note that therefore Eq. (13' 
is exactly equivalent to the original expression Eq. 
Reformionizing as described in Appendix ^ gives 

H,, = + r*V^(0)) b, (37) 



where we again use the Majorana representation for the 
impurity spin. The Hamiltonian Eq. (|3^) is very similar 
to the resonant level model that occurs in the solution 
of the two-channel isotropic Kondo model at the Emery- 
Kivelson line (compare Eq. ( p3[ ) with Jy — 0). The im- 
purity correction to the ground state energy can again 
be obtained from the Green's function. The propagator 
for the impurity has exactly the form Eq. ( p^ with F/j 
replaced by 



n. 



(38) 



The second equation can be verified by explicit calcu- 
lation of the sum of the two resonance energies given 
in Eq. (p5|). The impurity correction can be found from 
Eq. ( p6| ) in the channel-symmetric limit. Furthermore the 
density of states of the impurity was defined in Eq. ([2^). 
Combining these expressions we find for the impurity en- 
ergy 



2tt V 2r2 



TT V r 



(39) 



In the second equation we have used F <C Ec- It is now 
clear that the use of the two-channel anisotropic model 
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is justified whenever max{rQ,rb} « Fq + F;,. This con- 
dition is certainly met when N is very close either to an 
integer or to a half-integer value. Furthermore the range 
of A^- values for which the equation is approximately true 
will be larger for stronger asymmetry between Fa and F^ 
i.e. for a weak asymmetry between the reflection ampli- 
tudes jril and |r2|. 

In the two cases of electrons with spin in zero magnetic 
field or with a weak magnetic field we have max (Fa) <C 
max(Ff,) and the channel- anisotropic Kondo model is 
more convenient (See Fig. 2): In particular, the asymme- 
try simply produces a new energy scale obviously affect- 
ing the properties of the system close to the degeneracy 
point N = 1/2. This also allows to make explicit links 
with the small transmission limit (See Sec. VII). 



B. Applications 

We will now discuss a case for which the isotropic 
model is particularly well suited. This is the case of 
strong asymmetry between the reflection amplitudes of 
the two transmitted channels. We assume that one of 
the two channels is very close to perfect transmission, 
e.g. |ri| and |ri| <C \r2\ «C 1. Such a situation 
can be reached in a strong magnetic field where the elec- 
trons are essentially spin polarized. The number of open 
channels and their reflection amplitudes can then be ad- 
justed changing the voltage applied to the gates used to 
define the QPC. The wave-functions for the electrons in 
the two channels are '^a=i = ^"(x) ^^"^ ^q=2 — 
In the limit of interest here the energy Se with the help 
of Eqs. (H) and {Mj is found to be 



Se =. ^|r2p (1 + 2Acos(27rAf)) 
X ln((7|r2|V7r)(l + 2Acos(27r7V))) 



(40) 



where A = kil/1^21 is a small parameter. To leading 
order in A the correction to the capacitance is given by 

6C = +87£;c/3Vilk2|ln (|r2|^) cos{2ttN). (41) 

First we observe that at perfect transmission |ri| = 
there is dq ^signature of Coulomb blockade left as it is 
expectedO'B'tjEj. Furthermore it is interesting to notice 
that with one channel very close to perfect transmission 
the Coulomb blockade oscillations are strongly reminis- 
cent of the one-channel case, also discussed in Ref. ^. 



V. RENORMALIZATION GROUP 
FORMULATION 

We will now show that the energy scales Fa and Ff, 
can also be found from a renormalization group treat- 
ment. To do this we will use a similar argument as 
was put forward in Ref. n3 which in turn is based on 



Ref. g9[ We will first discuss the channel symmetric case 
kil — I ''2 1 (electrons with spin and zero magnetic field). 
The backscattering Hamiltonian for this case is given by 



y/ajEcv F cos(7r7V) cos(\/2^</>^ (0) ) . (42) 



The basic idea of the renormalization group treatment 
applied here is to calculate the partition function to sec- 
ond order in the small parameter |ri| and to choose |ri| 
as a function of the lattice step a in such a way that the 
partition function remains invariant under the transfor- 
mation a — > a' = a exp(Z). The correction to the partition 
function due to the backscattering perturbation is 



/3 

SZ^ - I dn dT2 {Hhs iTi)Hts (r2 )) 

Iq 



dTidT2 — ■ 

Vf\ti-T2\ 



(43) 



(27ra)2 



Using the definition of given in Eq. (|T^) it can be 
seen that the partition function 5Z is independent of 
the lattice step a and thus |ri| is invariant under rescal- 
ing. We now introduce a dimensionless parameter |fi| via 
1^1 1 = Jx/vf- With the help of this parameter we can 
rewrite the backscattering Hamiltonian in the standard 
form 



Hh. 



VF\ri\ 
2TTa 



cos(V27r0^(O)). 



(44) 



From this we can see that |fi| = A\ri\^J^Eca/ vp cos{ttN) 
has the meaning of an effective reflection amplitude. 
Since |fi| cx -v/akil it is clear that |fi| grows under the 
renormalization a — > a' = aexp(/). The renormalization 
flow equation for |fi| is given by 



(OI = ^k~i(OI 



(45) 



We can now integrate this equation from / = (a' = a) 
to I — Ic (a' = flc = aexp(Zc)), where Oc is the value 
of a' for which oc |fi| departs to strong coupling, 
that is |?^i(a')| ^ 1. Integrating we get |fi(ac)P = 
(ac/a)|fi(a)p = 1. The corresponding critical energy 
scale is defined through 



Ex,, 



vf/uc = VF\fi{a)\'^ /a = Jx"^ /{avp)- 



(46) 



Comparing Eq. ( p6| ) to Eq. ( |25| ) we see that E^^c is equal 
to Ff, up to a numerical factor. Further evaluation of 
Eq. (|46|) gives 



= 16|ri|27^cacos2(^iV). 



(47) 



Similar considerations can be used in the more gen- 
eral case where |ri| ^ |r2|. Then we have to use e.g. 
the backscattering Hamiltonian Eq. (|l6|) as a starting 
point. Expanding the partition function to the sec- 
ond order in this perturbation (the small parameters are 
\R\ = \ri\ + \r2\ and \5r\ = \r2\ - \ri\) we get 
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T 2 I T 2 

(27ra)2 



dTldT2 j — 

vf\ti 



^2 I 



(48) 



where the partition function is again independent of the 
lattice step a. Furthermore it is important to note that 
terms proportional to and Jy do not mix in the ex- 
pansion of the partition function up to second order. 

For a small channel anisotropy, either the coupling 
Jx or Jy grows under renormalization which gives rise 
to two independent energy scales. The critical energy 
for which {Jy) departs to strong couphng is E^.c = 
Jx'^/ivpa) (X Ec\R\^cos^{TrN) {Ey^c = Jy'^/ivpa) oc 
|(5rp-Ec sin^{7rA^)). Of course we have recovered here, 
up to a numerical factor the two energies ra(cx i?y,c) and 
r6(cx Ex,c) (see Eq. (H)). More generally, using Eq. (^) 
we find: 



d 
Jl 



(|i?(OP + IMOP)-(l^(OP + IMOP), (49) 



where we defined \R\ — Jx/vf and \Sr\ = Jy/vp- In the 
case of a strong channel anisotropy, we deduce that the 
effective coupling Jx^ + Jy^ flows off to strong coupling at 
the critical energy Ec = {Ex.c + Ey,c)- Again, we recover 
EcOcT = Th + Ta. For a small anisotropy, of course this 
reduces to E^ oc maxjFa, Ff,}. 



VI. REALIZATION IN A TWO-CONTACT 
SET-UP 

Let us now briefly emphasize that the channel- 
anisotropic two-channel Kondo model of Sec. Ill is also 
well suited to discuss the conductance behavior of a two- 
contact set-up a la Furusaki-Matveev. In this setup 
which is a natural extension of the geometry shown in 
Fig.|I| a single dot is coupled to two different electron 
reservoirs via two point contacts. A back-gate or side- 
gate is used to vary the Coulomb energy. The device is 
illustrated in Ref. |l^. 

Let us start with a strong magnetic field such that 
the electrons of each reservoir are fully polarized. We 
obtain an effective two-channel model where a = 1(2) 
now denotes the electrons of the left (right) QPC. Again, 
we are mainly interested to the case of a small channel 
anisotropy where reflection amplitudes at each contact, 
namely |ri| and \r2\, are close to each other: |ri| < |r2| <C 
1 . It is convenient to associate the centers of the two con- 
strictions with the point a: = 0, meaning that electrons in 
the dot are described by '^i^l{r) at a; > and by '^2,l{r.) 
at a; < 0li3. We introduce the variables (jjt and (jjc via 



(50) 



and similarly for 9i and ^2- Typically, here (/)c is the 
bosonic variable associated to the charge in the dot, while 
(j)t is the one related to the current passing through the 
dot. The main difference with the one-contact set-up is 



that here the charge on the dot rather reads (note the 
minus sign)E3 



Q 



= (0l(O)-02(O)). 



(51) 



It is straightforward to show that the backscattering 
Hamiltonian is of the form of Eq. (16) or Eq. (18) with 
(ps replaced by explicitely 

Hbs = ^ cos (^V2^M0)) ^^ + ^ sin (V2^</)t(0)) Sy. 

(52) 

The Kondo parameters are given in Eq. ( |l9[ ) where 
\Sr\ ^ \R\ <C 1, and the kinetic energy for the sym- 
metric charge mode HKin{<Pt,St) is still of the form of 
Eq. (||). Again, we insist on the fact that this mapping 
is more intuitive than the one of Eq. ( [l7| ) because the 
anisotropy of the reflection coefficients is directly related 
to the anisotropy between coupling parameters in the 
Kondo model. In consequence, we have two independent 
energy scales Fo and F;,. The behavior of the conduc- 
tance close to the degeneracy point A'' = 1/2 becomes 
really transparent. 

In the absence of any backscattering, the system is 
equivalent to two resistances 2TTh/e'^ (those of the two 
QPCs) connected in series. Therefore the conductance is 
Go = e'^ / {AttTi) . This is in fact still the case for finite bar- 
riers p, the symmetric case |ri| = \r2\ at the degeneracy 
pointE2l, where both Jx and Jy are zero for A'' = 1/2. As 
was found in Ref. |l^ the conductance is still at its reso- 
nant value Go and the tails of the peak at A^ = 1/2 are 
not Lorentzian: G (x (A^ — 1/2)^^ due to the Jx coupling. 
When \Sr\ is finite, from Eq. (52) we immediately see 
that the asymmetry between channels engenders a finite 
backscattering process at A^ = 1/2 (again, Jy cx \Sr\): 



J y 



bs 



3in(V2^0t(O)) Sy, 



(53) 



This produces the same physics as a non-magnetic im- 
purity in a Lultinger liquid with the Luttinger exponent 
being g — 1/2l3. The problem becomes exactly solvable 
using the refermionization procedure of Appendix A. In 
particular, tbe current through the device takes the re- 
quired formed: 



evFil^HO)tp{0). 



(54) 



For temperatures T ^ Ey^c and A^ close to 1/2, the ef- 
fective potential scattering Jy (or the asymnietxx param- 
eter) diverges and the conductance then obeyfj: 



G(T;A^«l/2) = Go 



2/3-2 



T 



ErN"^- 



(55) 
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where the second equation is true for g = 1/2. The en- 
ergy Ey^c was defined in the last paragraph of the previ- 
ous section. For non-symmetric barriers the conductance 
peak height (like the capacitance peak height before) be- 
comes strongly dependent on the finite asymmetry be- 
tween reflection amplitudes at the QPCs (Fig. 4). This 
naturally reflects the fact that the on-resonance behav- 
ior of the system for TV = 1/2 is very different if the two 
barriers are not identical. 

A simple explanation can be given using the channel- 
anisotropic two-channel Kondo model. On resonance, 
the small asymmetry inevitably grows under renormal- 
ization. At low temperatures T <C m.ax.{Ta,Yb\, the 
conductance behavior for all values of N is given by 



G{T-N)^Gc 



T 



This can also be rewritten a 

'T 



< Go. 



G{T-N)^{- 



(56) 



(57) 



where (See Eq. (38)) 

Fcx (|ri|2 + |r2p + 2|ri||r2|cos(27rA^)), (58) 

even though this formula seems a little bit less intuitive. 
In particular, the N-dependence of the conductance for 
iV w 1/2 is less apparent. It is worth noting that the 
quadratic temperature dependence is a universal prop- 
erty of inelastic co-tunneling. But, close to = 1/2 this 
can also be interpreted as a nice manifestation of the 
restoration of the Fermi-liquid behavior due to the finite 



0.1 
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FIG. 4. Conductance in the two-contact set-up (in strong 
magnetic field) as a function of A'" for diflFerent values of the 
(small) anisotropy parameter \Sr\ = \r2\ — \t\\. The parame- 
ter \R\ = In I -I- jr2| is here set to \R\ = 7/10. The diflFerent 
curves have been calculated for T/Ec = 1/50 and have been 
normalized to Go. In order of decreasing conductance peak 
height the diflFerent curves here correspond to \Sr\ = 0.2 (top), 
\5r\ = 0.25, \5r\ = 0.3, \Sr\ = 0.35 and \5r\ = 0.4 (bottom). 



asymmetry between channels (like for the capacitance 
problem). 

Unfortunately, the present approach does not allow a 
complete solution for the case of fermions with spin in 
the two-contact set-upEZl. At zero magnetic field and for 
completely symmetric barriers at the two QPCs, there 
exists a mapping onto a four-channel Kondo model for 
the Hamiltonian (which could help to compute thermo- 
dynamical properties) , but it is very difficult to rigorously 
compute current-current correlation functions in the new 
basis and then to extract the conductance behavior; The 
current operator indeed exhibits an unusual form as in- 
dicated in Appendix B of Ref. [l^. Therefore, the height 
of the conductance peak cannot be calculated in a rig- 
orous manner even for completely symmetric barriers at 
zero temperature. Furthermore, the mapping onto the 
Kondo problem is only valid at (close to) N = 1/2. To 
extend the (four-channel) Kondo mapping to the case of 
(weakly) asymmetric barriers remains a challenging task. 



VII. SMALL TRANSMISSION LIMIT 



Having so far concentrated on the limit where the 
quantum dot is strongly coupled to a reservoir through 
a highly transmissive quantum point contact, we will in 
this section consider the limit of weak coupling. The 
similarity between the Coulomb blockade problem in this 
limit and. a Kondo model was noticed by Glazman and 
Matveevul. An explicit mapping of the Coulomb blockade 
Hamiltonian on a Kondo model was used by Matveev to 
calculate the charge and the capacitance of the quantum 
dot in the weak transmission limitEl. Note that recently a 
noncrossing approximation (NCA) has been, generalized 
to this type of multichannel Kondo modelsEil. 

We will here rederive Matveev's mapping and discuss 
the straight-forward extension of his model to the channel 
anisotropic case. Again our discussion will be restricted 
to the case of two transport channels through the point 
contact. The anisotropy between the transmission ampli- 
tudes for the two transport channels in the QPC will still 
give rise to a mapping on a channel-anisotropic Kondo 
model; But here, the system flows off to the usual spin- 
isotropic fixed point, i.e., Ja,i_ — Ja,z ^ 0. We will see 
below that the problem of calculating the average charge 
(Q) on the dot is now equivalent to the problem of find- 
ing the average of the z-component {Sz) of the impurity 
spin in the Kondo model. In addition the capacitance 
of the dot is the same as the magnetic susceptibility x 
of the impurity. Once this equivalence is established the 
problem is solved, since these quantities {{Sz), x) can be 
found in the literature. The channel- and spin isotropic 
multichannel Kondo problem was solved exactly in this 
limit in Refs. while the impurity-susceptibility for 

the channel-anisotropic (but spin-isotropic) two-channel 
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Kondo problem can be simply extracted from Ref. |3J. 
Note in passing that the channel-anisotropic case was also 
solved (exactly) using Bethe Ansatz in Ref. There, the 
Wilson ratio is computed even in the case of a channel 
anisotropy. This model has also been investigated using 
conformal field theory and numerical renormalization- 
group calculations in Ref. |36| . 

We will now proceed to writing down the model for 
our system in the small transmission limit. Instead of 
formulating it in momentum space as it was done by 
Matveev we will here use a real space formulation which 
allows us to stress the analogy with the corresponding 
model in the strong transmission limit (see Sec. ||). In 
Appendix we will also address the question why the 
mapping on the Kondo model can not be derived using 
bosonization as it was done in Sec. In the previ- 
ous section we have treated the backscattering as a small 
perturbation to an otherwise perfectly transparent QPC. 
In this section the perturbation is a tunneling Hamilto- 
nian which couples two a priori independent systems (the 
2DEG and the dot). A smooth transition can be made 
from the strongly to the weakly coupled limit by con- 
tinuously increasing the auxiliary gate voltage to pinch 
off the QPC. In this transition a perfectly transmissive 
one-dimensional channel n will be cut into two weakly 
coupled halves. In the vicinity of the center of the QPC 
electron motion is still quasi one-dimensional. Electronic 
wave-functions to the left of the center of the QPC (in 
the reservoir) will be denoted q, wave- functions to the 
right of the center (in the QD) are ^'^ i- Again a is the 
spin index, n the channel index due to lateral confinement 
and here a = 0, 1 indicates the location of the electron. 
Hopping between the two sides of the QPC constitutes a 
small perturbation. To model this perturbation we use 
the Hamiltonian 



Hj 



E 

cr=l,2 



(|i,|<i(0)*,,o(0) 



h.c. 



(59) 



Note that the kinetic energies Hfej„ \ a.a}^ and 

^kin for electrons in the 2DEG and the 



dot have the form of Eq. (|22|). The boundaries for the 
integration along the x-axis are — oo,0 in i?{!j„ and 0,oo 
in Hh^. The Coulomb interaction can be modeled as in 
Sec. Q), Eq. (^ the charge on the dot being 



Q 



dx 



(^'t,i(x)*,,i(a;)-po 



(60) 



The equilibrium charge density po is chosen in such a 
way that the total charge Q on the dot is zero when no 
voltage is applied to the gate (Vg — 0). 

We now want to concentrate on the point N = 1/2 
where the states with Q = and Q ~ 1 are energy 
degenerate (See Eq. (||)), and therefore charge fluctu- 
ations are large. We introduce the small parameter 
U — e/{2Cgd) — Vg cx (A^ — 1/2) to measure deviations 



of N from the degeneracy point. In terms of U the elec- 
trostatic energies of the two states Q = and Q = 1 are 
Eq = Q and Ei = eU respectively. 

If the two conditions \U\ <^ e/Cgd and kT <C Ec are 
met only the two states with (5 = and Q — 1 are acces- 
sible and higher energy states can be removed from our 
theory introducing the projection operators Pq and Pi. 
Here Pq and Pi are projecting on the states Q — Q and 
Q — 1 respectively. The effective Hamiltonian for this 
truncated system is 

Heff = {Hl^ + Hl^) (Po + Pi) + eUPi (61) 
+ E (|i.|<i*..oPo + |<a|<o*-.iA), 

cr=l,2 

where the operators in the tunneling part are evaluated 
at a: = (see Eq. ^). From this equation we see that 
dH/dU — ePi and therefore {Q) ~ dEo/dU where Eq is 
the ground-state energy of the Hamiltonian Eq. (|6l|). It 
can be shown that the Hamiltonian Eq. ( |6l| ) is equivalent 
to the Hamiltonian 



= E ^fe- + 2^^- (62) 
Q=oa 

+ ^ E {J..± {4(0)3- + s-{0)S+)) 

(7=1,2 

which has the form of a standard Kondo Hamiltonian. 
The equivalence between the Hamiltonians Eq. (0) and 
Eq. ( |6^ ) was demonstrated by Matveev in Ref.^ For 
completeness we rederive this relation in Appendix ^ 
taking channel anisotropy into account explicitely. In 
Eq. (^) we have introduced the transversal Kondo cou- 
pling parameters Ja.i_ — '^\t(7\ and an effective external 
magnetic field h — eU/2 oc {N — 1/2) applied along 
the z axis. Note that Eq. ( |6^ ) corresponds to the limit 
Ji.z — 'h.z = of Eq. (|l|). We have dropped the con- 
stant eU/2 from the Hamiltonian in Eq. (|62|). Therefore 
the gro und state energy Ek^o of the Kondo Hamiltonian 
Eq. ( p^ ) is related to the ground state energy Eq of the ef- 
fective Hamiltonian Eq. (61) through Ek^ + eU /2 = Eq. 
Since 2{Sz) = dEKfl/dh and (Q) — dEQ/dU we are led 
to the obvious identification 



1 



{Sz 



(63) 



Combining Eq. (g) with U — e/{2Cgd) — Vq we find for 
the correction to the capacitance C = —d{Q) / dU (x x — 
d{Sz)/dh, where x is the impurity susceptibility (Since h 
only acts on Sz, X is also equivalent to the local magnetic 
susceptibility xi)- 

So we have now shown that a calculation of the capac- 
itance of the quantum dot coupled to a 2DEG through 
a QPC with two transport channels with different trans- 
mission amplitudes is equivalent to a calculation of the 
impurity susceptibility in the channel-anisotropic two- 
channel Kondo model. It is known that both J^^z = 
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and J„^± <C 1 grow under renormalization and that for 
small enough energies Ja,z = Jcr,i^ = Jo- oc \t^\ even if we 
start from a spin-anisotropic Kondo model. The resulting 
spin-isotropic Kondo model was solved in Refs. |3^,^. In 
the zero temperature limit the impurity susceptibilities 
for the one-channel case and the two-channel isotropic 
limit are 



est. J2 = 0, 
\ogh Ji = J2. 



(64) 



Note that the susceptibility for the two-channel case di- 
verges in the zero magnetic field limit. The channel- 
anisotropic (but spin-isotropic) Kondo. model was dis- 
cussed e.g^by Coleman and SchofielcCJ (also by Andrei 
and JerezEa and by Affleck et alE3) who found for the 
susceptibility at zero magnetic field (h — 0) 



X (X 1 - 



\og{iy) 



(65) 



where the anisotropy parameter v oc \St\/T with \6t\ oc 
I Ji — J2I and T oc (Ji -I- J2), is z/ = 1 in the one-channel 
limit (J2 = 0) and v = in the two-channel isotropic 
limit. Eq. ( |65| ) thus reproduces correctly the main char- 
acteristics of the result Eq. (|6^ ) in these two limiting 
situations. From the results for the susceptibility we im- 
mediately obtain the capacitance 



6C< 



est. J2 = 0, 

log(7V-l/2) Ji = J2, 
\og\St\ Ji^J2,N 



1/2. 



(66) 



The results for the one-channel aad the two channel 
isotropic cases are due to MatveevB. The main obser- 
vation we want to make here is that a small channel 
anisotropy cuts off the logarithmic divergence exactly in 
the same way as in the case discussed in the previous sec- 
tions where the reflection amplitudes in the QPC could 
be treated as small parameters. 



VIII. CONCLUSIONS 

We have applied the channel- (and spin-) anisotropic 
two-channel Kondo model to study Coulomb blockade os- 
cillations in the capacitance of a quantum dot. Our main 
interest has been to investigate the effect of an asymme- 
try between the reflection (or transmission) amplitudes 
of different open channels in the QEp connecting the dot 
to a reservoir. Following MatveevHcl we have studied the 
two exactly soluble limits of very weak (|ti|, 1^2! ^ 1) and 
very strong coupling (|ri|, \r2 \ ^ 1). 

A summary of the results for the capacitance in differ- 
ent limits is given in Table 1. 

In both limits a mapping of the original problem onto 
a Kondo model is possible. Remember that these re- 
sults concern the limit of very low-temperature where 



quantum fluctuations are prominent. At quite high tem- 
perature, thfj-^Cojrdo physics gets destroyed by thermal 
fluctuationifiy^N. 

For weak backscattering at the point contajet 
(ki|,k2| ^ 1) the original problem can be mappedtil 
on a chaanel-anisotropic Kondo model at the Emery- 
Kivelsor£3 line (Ji,z = J2,z = 27rwi?). For this partic- 
ular value of the coupling constant the Kondo model is 
exactly solvable. The anisotropy of the reflection am- 
plitudes for different channels is directly reflected in the 
channel-anisotropy of the Kondo model. In fact we found 
Ji,± oc |ri| and J2^± oc |r2|, where Ji^± and J2,± are the 
coupling constants for the two channels in the Kondo 
model. The mapping allowed us to calculate the shift of 
the ground-state energy of our Hamiltonian due to the 
backscattering at the point contact and in turn to find 
the capacitance of the quantum dot. While the capaci- 
tance is logarithmically divergent for values of the gate 
voltage close io = n -I- 1/2 (ri is an integer) in the 
isotropic limitlj (|ri| — \r2\), this divergence is cut off by 
a small anisotrop«-.between the reflection amplitudes of 
different channelsO. 

This can be interpreted as a manifestation of the 
restoration of the Fermi-liquid behavior close to the de- 
generacy points N = n -\- 1/2 due to the asymmetry be- 
tween channels. 

Note that a similar conclusion can be reached investi- 
gating the conductance behavior in a two-contact set-up 
in a strong magnetic field where |ri| and \r2\ denote the 
backscattering amplitudes at the two different point con- 
tacts. The on-resonance behavior (G = Gq = e^/2/i) for 
A^ = 1/2 is reduced for a small anisotropy between |ri| 
and \r2\. 

There are two intrinsic energy scales (F;, oc (|»'i|-|-|r2|)^ 
and Fa oc (|ri| — |r2|)^) in the channel-anisotropic Kondo 
model. Each can be interpreted as the resonance en- 
ergy related to coupling half an impurity to the conduc- 
tion electrons. In the channel- isotropic case only half the 
impurity spin is screened by the conduction electrons. 
Coupling back the second half of the impurity to the 
conduction electrons leads to the emergence of a second 
energy scale Fq. It is this new energy scale that en- 
ters the expression for the capacitance and cuts off the 
divergence at A^ = 1/2. Unfortunately there is no di- 
rect correspondence between the magnetic susceptibility 
in the Kondo model and the capacitance of the quantum 
dot. Such an equivalence exists only in the limit of small 
transmission. It is true however that the behavior of the 
capacitance (found from the Kondo model at the Emery- 
Kivelson line) is reminiscent of the impurity susceptibility 
X = d'^{6e)/dh^ away from the Emery-Kivelson line or of 
the local magnetic susceptibility xi = d{Sz)/dh at the 
Emery-Kivelson line. 

In this paper we have extended our previous workEl 
deriving the mapping on the anisotropic Kondo model in 
a pedagogical way and carefully discussing its limits of 
validity. We have then given an alternative way for calcu- 
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lating the charge on the quantum dot using a mappingEzI 
on a channel-isotropic two-channel Kondo model. In this 
approach the coupling constant is a complex parameter 
depending on both |ri| and \r2\- While the latter ap- 
proach has the advantage of being exact, it seems less in- 
tuitive since the anisotropy of the reflection coefficients is 
not directly reflected as an anisotropy between coupling 
parameters in the Kondo model. We have in addition 
used a simple scaling argument to recover the intrinsic 
energy scales Ta and Yb that occur in the expression of 
the capacitance. 

We relied on purely mathematical arguments to show 
the equivalence between the Coulomb blockade problem 
and the Kondo Hamiltonian in the strong tunneling limit 
(ki|,k2| ^ !)• In the opposite hmit |ti|,|t2| ^ 1 the 
similarity of these two problems can be understood us- 
ing a comparably simple physical argument (see Ref. ^ 
and Appendix The main observation is that at low 
enough temperatures T <^ Ec and for voltages close to 
A'' = n -|- 1/2 only two charge states on the quantum dot 
are energetically accessible (e.g. Q = 0, Q = 1). The 
charge state of the dot is then interpreted as a pseudo- 
spin-1/2 degree of freedom which corresponds to the im- 
purity spin in the Kondo model. The real spin of the 
conduction electrons in the Kondo model is replaced by 
an index a indicating the location of an electron (distin- 
guishing between electrons to the left and to the right of 
the QPC) in the Coulomb blockade problem. 

A first order process in the Kondo problem that flips 
the impurity spin from up to down and the spin of a 
conduction electron from down to up is then equiva- 
lent to a tunneling process that takes an electron from 
the left to the right and changes the charge on the dot 
from Q = to Q = 1. The Kondo coupling parame- 
ters simply are Ji(2),± |ii(2)| and Ji(2).z — 0. As it 
is clear from the derivation of the mapping there exists 
an equivalence between the charge on the dot and the 
impurity spin, namely (Q) cx (Sz)- It was furthermore 
shown that the capacitance is basically the same (up to 
some constant) as the impurity susceptibility. Since the 
model with Ji[2),± ^ l*i(2)l and Ji(2),z = flows to 
the usual spin-isotropic fixed point we have been able 
to use the known result for the impurity susceptibility in 
the channel-anisotpopic two-channel Kondo model with 
Ja,i_ = Ja.z S> OEJ, to discuss the effect of channel- 
anisotropy on the capacitance. Exactly as in the limit 
of strong transmissioa the capacitance diverges in the 
channel-isotropic cascEl, but the divergence is cut by the 
anisotropy. Note that no immediate connection can be 
made with the Kondo model at the Emery-Kivelson line 
("^1(2), 2 = '2ttvf) which we obtained when treating the 
small reflection limit. 

Let us finally remind that all the results on the behav- 
ior of the capacitance of the dot close to A^ = 1/2 have 
been summarized in Table 1. 
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APPENDIX A: REFERMIONIZATION 

In this Appendix we want to elaborate on the unusual 
refermionization procedure which we used in Sec. (III). 
The backscattering part of the bosonic Hamiltonian to 
be refermionized is given in Eq. dH). The kinetic energy 
is of the form 



Hkv 



dx 



(Al) 



It will turn out to be convenient to use the field tt^ (x) in- 
stead of 9s{x) for the moment. The two fields are related 
through TTsix) = dxOs{x). The commutation relations for 
the fields (j)s{x) and tTs{x) are [4>s{x),iTs{y)] = i^{x — y), 
[(j)s{x),(j)s{y)\ = and ['Ks{x),'Ks{y)\ = 0. The basic idea 
of the refermionization procedure is to introduce an op- 
erator ^{x) such that cos(\/27r(/)s(0)) (x ?A(0) -t-?A^(0) and 
sin(V27r0s(O)) oc V'(O) ~ V'^(O). It is clear that for such 
an operator 



V'(O) 



exp(iv27r(/)s(0)). 



(A2) 



In addition the operator must obey the usual fermionic 
anti-commutation relations. Using the relation e"^e^ = 
e^e^e^^'^^ it can be seen that the obvious choice for t/i, 
namely — (27ra)~"'^/^ exp(zV27r0s(x)) does not obey 
anti-commutation relations. To construct a fermionic op- 
erator we introduce the auxiliary fields 



(t)±{x) = {(t>s{x) ± 0s (-a;)) , 
7i"±(a;) ^ {n^ix) ± TTsi-x)) . 
In terms of these new fields the kinetic energy is 



(A3) 



H 



Ki 



VF J2 

a=± 



dx 



{dx4>af + {dxOaf 



(A4) 



We thus arrive at a theory which is confined to positive 
values of x. The advantage of this restriction becomes 
clear when we introduce the two additional right-going 
and left-going fields 



^ R..±{x) = (t)±{x) - I dy Ti±{y) 



(A5) 



^L,± (a;) = (t)± {x) + dyTT± (y)) 



At x = we have $fl,+(0) = $l,+(0) = V2(j)s{0) and 
also ^R-{0) = (0) = which makes these fields 
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candidates for the construction of our new fermions. The 
backscattering Hamiltonian can be expressed through the 
fields ^R^+{x) and $l_+(x) only. There are in fact many 
ways in which this can be done. However, we will soon 
get rid of this ambiguity. The fields ^r{l),-{x) occur 
only in the kinetic part of the Hamiltonian and are thus 
of little interest to us. Later on we will need the commu- 
tation relations 



operators we have defined above the kinetic energy fi- 
nally takes the simple form given in Eq. ( ^2|) while the 
backscattering part of the Hamiltonian is 



Hbs — 



Jx 



(A14) 



=+isgn(a;-y), 
[^L,+{x),^L,+ iy)] ^ -isgn{x- y). 



(A6) 



For completeness we also give the correlation functions 
Gr(l) = ($H(L),+ (a;)*fl(L),+ (0) - $i?(L),+(0)2) which 
take the standard form 



Gr(l) = -In ( ° . 



(A7) 



The plus sign belongs to the label R while the minus sign 
belongs to L. The kinetic energy in terms of these new 
fields takes the form 



HKin — 



VF 



2 .ti 



dx 



(A8) 



We now drop the $fl(i)^_(a;) part of the kinetic energy 
since it is not coupled to the backscattering term. To 
refermionize the ri^l) ^+{x) part we introduce the oper- 
ators 



^r{x) 
iI^l{x) 




iy7F$H, + (a;) 



(A9) 
(AlO) 



Using the commutation relations for the bosonic fields we 
can verify that the fields '>jjR^±{x) and ipL,±ix) really are 
fermions and obey 



(All) 



where p = R,L. No te, t hat for a; = we have ipL.+ {0) — 
i/'i?,-i-(0)) (see Eqs. (A5) and ([A9|)) and there thus seems 
to be more than one way to refermionize the backscat- 
tering Hamiltonian Eq. (jl^). To lift this ambiguity we 
reextend our theory on the full a;-axis via the definition 



ip{x) = V 



'4)r{x) X > 0, 
V-'l(— a;) X < 0. 



(A12) 



In the above definition of the fcrmion 'ip{x) (Eq. ( [A12| )) 
we have introduced an additional phase factor 



r = exp{iTTdU) = l-2dU 



(A13) 



to ensure that ip{x) anticommutes with the spin operators 
Sx and Sy written in terms of Majorana fermions d and 
d'^ (cf. Eq. (|l|)). The second equality in Eq. ( |A13D holds 
because d^d = 0, 1 at zero temperature. In the fermionic 



Using Eqs. (A13) and (21) together with the commuta- 
tion relations for the operators d and d^ we can show that 
VSx — —iSy and VSy — —iSx- With these relations we 
recover Eq. (H). 



APPENDIX B: GREEN'S FUNCTIONS 

Although this is rather standard material (see e.g. 
Rcf. 11^) we believe that it is useful to give a short 
derivation of the Green's functions Eq. The Fourier 
transforms of the impurity Green's functions Gair) = 
-{Tra{T)a{0)) and Gfe(r) = - {Trb{T)b{0)) can conve- 
niently be found from the equations of motion. We will 
here only derive the correlation function for a, the cor- 
relator for b can be found along the same lines. The 
Hamiltonian H^j^ = HKin + ^^6s of our system is given 
in Eqs^^ ^ and (H). The fields ij^ and d = a + ib (see 
Eq. (|2l|)) obey standard fermionic anti-commutation re- 
lations. Introducing the Majorana components of the 
field through zi{x,t) — {^\x,t) +tp{x,T)) /V2 and 
Z2(x,t) = {iP{x^t) — ij)'^ {x^rXyj we see already 
from the Hamiltonian in Eq. ( psj) that a couples only to 
Z2- We introduce the additional propagator Gz2a{x,T) — 
-~{TT-Z2{x,T)a{0)). The equations of motion for the two 
coupled correlators Gair) and 6*220(2;, r) are 

drGair) = -5{t) + iJ^Ca{Q, t), (Bl) 
V27ra 

drGz2a{x,T) = +iv pd^G ^^ai^ , t) - i-^^==b{x)G aij) . 

To solve these equations it is best to go to Fourier space 
making use of the relations 



G.,a(x, ^) = ^ E / "-+'f-G,,Jp, c^„), (B2) 

where the sum is over the fermionic Matsubara frequen- 
cies ojn = (2n-|- 1)tt/ p. To calculate the correlator Ga(T) 
we only need to understand the local physics in x = 0. 
The local equations of motion in Fourier space are 



zw„Ga(w„) = 1 - i- 



J n 



G^^aM = i^^G^°\iUn)Ga{uJn)- 

vzTTa 



(B3) 
(B4) 
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To alleviate the notation we have introduced G^^\uJn) = 
— i sgn(w„)/2vi?. To find G'-*'-'(w„) we Fourier transform 
the free electron propagator G'^^'> (p, aj„) — {iujn — vfp)^ 
with regard to p and take the limit x — > 0. Furthermore 
we have defined G^^a(j£n) = G„a{x = O.ujn)- 



Substituting Eq. (B4) into Eq. (B3) we can solve for 
Ga{ijJn) and obtain 



Ga(Wn) 



1 



iuJn + iTaSgn{ljJn) ' 



(B5) 



After an analytic continuation icij„ — > a> + we recover 
Eq. (H). 



APPENDIX C: MAPPING TO THE KONDO 
MODEL IN THE SMALL TRANSMISSION LIMIT 

In this Appendix first we want to fill in the gaps be- 



tween Eq. (|6l|) and the Kondo Hamiltonian Eq. (f 



Let 



us consider the tunneling part of the effective Hamilto- 
nian. The first term takes an electron from the 2DEG 
and transfers it to the QD, the projection operator Po 
makes sure that the charge on the dot is Q = before 
the tunneling takes place. The second term takes an elec- 
tron from the dot to the lead. In our truncated system 
this process is allowed only when the charge on the dot 
is Q = 1- This restriction is implemented through the 
operator Pi. The main goal of the following manipula- 
tions will be to show the equivalence of the Hamiltonian 
Eq. ( |6l| ) to a Kondo Hamiltonian. 

To explicitly account for the charge on the dot we make 
the replacement |$) \^)\Q). Here 1$) is any state of 
our system with charge Q on the dot. The values of Q are 
limited to Q 0, 1 and the states \Q) = |0) and \Q) — |1) 
can be considered as the basis of a two-dimensional vector 
space. However, the product is no tensor product 

since the charge of the dot is of course not independent 
of the system's state. The state \Q) should rather be 
considered as an auxiliary label to |<I>). In addition to 
introducing the label \Q) we make the replacement 



(Cl) 



in Eq. (61). Here 5*+ and S are pseudo-spin ladder 



operators acting only on the charge part | Q) . 

Since S+\Q = I) = and S'-jQ = 0) =0 these op- 
erators ensure in the same way as the projection opera- 
tors Po and Pi that only transitions between states with 
(5 = and Q = 1 take place. In addition the charge 
Q on the dot is adjusted whenever a tunneling process 
takes place since S+\Q) = |1) and S'-|l) = |0). We would 
like to emphasize again that only the combinations of 
pseudo-spin ladder operators and hopping operators in- 
troduced above are meaningful since \Q) and |$) are not 
independent. 



To get rid of the remaining projection operators in 
Eq. (|ll) we rewrite {eUPi) as [eU {Pq + Pi) / 2 + eU {Pi - 
Po)/2]. We observe that 



(Pi±Po)|0) =±|0) 
(Pi±Po)|l) =+|l) 



(C2) 



This leads us to identify (Pi -Pq) with 2S^ and (Pi -f-Po) 
with the identity operator on the space spanned by |0) 
and |1). Here we used that for the z-component Sz of the 
pseudo-spin we have S'^|l) = |l)/2 and 5*^10) = -|0>/2. 
Gathering all terms we can rewrite the effective Hamil- 
tonian Eq. ( |6l| ) as 



H, 



{P^kir 

E 

(T = 1.2 



f Hl,^) + eU{2Sz + l)/2 (C3) 

t,\-^l,^^^0S+ + \t„\^l,^^aAS-) . (C4) 



We now introduce an additional pseudo-spin operator 
sj(x) via 



4 = o*.4 



(C5) 



where the matrices a ~ a^i: i(Jy are standard combina- 
tions of Pauli matrices. It is important to note that these 
pseudo-spin operators again have nothing to do with the 
true spin of the electrons but are related to the location 
of an electron (a = for an electron in the 2DEG, a — 1 
for an electron in the QD). Int rod ucing the definition 
of these pseudo-spins into Eq. (C3) finally leads us to 
Eq.®. 

Finally, we want to show that in this limit the 
bosonization approach does not allow us to precisely 
build a pseudo-spin operator describing the dot from 
the original tunnel Hamiltonian. For simplicity, we re- 
strict the discussion to the case of spin- less fermions, i.e., 
we ignore the spin index cr. The main problem we en- 
counter is that at a; = we have open boundaries, im- 
plying that *i(0) = *q(0) = 0. Introducing right (R) 
and left (L) movers as in Eq. (|^), this is equivalent to 
write, e.g. for the dot, 4'i,_r(0) -I- 5'i,l(0) = 0. This has 
the effect to pin the (charge) fields 0i and (po a.t x = 0: 
(/)i(0) = 0o(O) = V^/2. Therefore, this provides us 



(G6) 



From the form of the tunnel term (which can be rewrit- 
ten either with 'i'l.n or with ^^i^l) we would be tempted 
tOpj^xplicitly build the pseudo-spin operator in the dot. 



(C7) 



(0) - 1/2. 
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However, due to the open boundary condition at a; = 0, 
the fermion operator ^'i^jjjx at a; = now only depends 
on the superfluid phase Then, 5"+ would commute 

with Sz and then S would not be a quantum spin object. 
The only way to proceed in order to recover the (correct) 
Kondo mapping is to introduce the extra label \Q). 
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Nbr of channels 


Reflection or transmission coeflicient 


Capacitance at N=l/2 


Ref. 


1 


In 1 < 1 (\ro\ —f 1] 

I'll w z\ ^ } 


const. 


Ref. 


9 




2 


In 1 = Irol < 1 

1 ± 1 1 Z 1 


~ln\N - 1/2| 


Ref. 


9 




2 


\ri\ + \r2\ < 1 


-ln(|r2| - |ri|) 


Ref. 14, Sees. [II V 


1 


« 1, (|t2| ^0) 


const. 


Ref. 


5 




2 


\ti\ = 1^21 < 1 


-ln|iV- 1/2| 


Ref. 


5 




2 


|ii|^|i2|«l 


-ln(|<i|-|t2|) 


Sec. \ 


a 


I 



Table 1: We have here listed all the results for the capacitance of the quantum dot at N ^ 1/2 and the respective 
references. The divergence occuring at N ^ 1/2 in the channel-isotropic two-channel limit is cut off by a channel 
anisotropy for both weak and strong reflection at the QPC. 
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